**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 

Radiatively Driven Winds of OB Stars — from Micro to 
Macro 

Jifi Krticka 

Ustav teoreticke fyziky a astrofyziky PfF MU, Kotldfskd 2, 
CZ-611 37 Brno, Czech Republic, email: krticka@physics.muni.cz 

Jifi Kubat 

Astronomicky ustav AV CR, CZ-251 65 Ondfejov, Czech Republic, 
email: kubat@sunstel.asu.cas.cz 

Abstract. We review basic physics of line-driven stellar winds of OB stars. 
We discuss elementary processes due to which stellar winds are accelerated on a 
microscopic level. We show how these microscopic processes may enable the out- 
flow and how they determine wind properties on a macroscopic level. We discuss 
shortcomings of present wind theories and future wind model improvements. 



1. Introduction 

Observations of many OB stars show that there exists an outflow of material 
from the stellar surface into the interstellar medium - the stellar wind. The the- 
oretical study of hot star winds started few years after the discovery that elec- 
tromagnetic radiation carries mo mentum tha t ca n be transferred to th e matter 
in the process of light scattering. iMilnel (1192.^ and l.Tohnsonl (Il925l.ll92fil) stud ied 
the possibility of the emission of high-speed atoms from stars. iMilnd (|1926h in 
a beautifully written paper showed the importance of the Doppler effect for the 
line radiative acceleration. However, for the next few decades the radiatively 
driven hot star winds did not attract much attention. Modern studies of hot 
stars' wi nds were stimulated main ly b y UV observations of hot stars. P ioneering 
works of iLucv fc Solomonl (|l97tf l and ICastor. Abbott fc KleU (|197,1 hereafter 
CAK) serve as a basis for present h ot star wind theory. 

Except a more general book bv lLamers Sz Cassinellil (|l999l ). lOwockil (I2001T) 
provid es an elegan t intro duction to the hot star wind physics, while lKudritzki &: Pulsl 
( 200(1 or lOwockil (|2004T t can be consulted for a more detailed review. 

In this review we intend to answer two basic questions: First, what are 
the processes responsible for the acceleration of the stellar wind of OB stars? 
Second, what are the theoretical predictions of the wind structure? 

2. Stellar wind of hot stars: micro- view 

According to the theoretical studies of hot star winds, they are accelerated due 
to the light scattering in lines of heavier elements, and to a lesser extent due 
to scattering on free electrons. Hydrogen and helium are mostly inefficient for 
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Figure 1. A typical wind volume with 1000 H ions, 100 He ions, 2 metallic 
ions, and corresponding number (about 1200) of electrons. 

wind driving. From the microscopic view of stellar wind in Fig. ^ it is clear that 
two basic processes are necessary to accelerate the bulk wind flow: 

1. process that transfers momentum from the radiation field to heavier ions 
and electrons, 

2. process that transfers momentum from heavier ions and electrons to the 
bulk flow (hydrogen and helium ions - mostly passive component). 

2.1. How to transfer the radiation momentum to heavier ions? 

To better understand transfer of the radiation momentum to heavier ions we 
shall generally discuss absorption and emission mechanisms. The photon fre- 
quencies before absorption and after emission are v and i>', respectively (note 
that the case v 1 = means pure absorption). 

In the general case v ^ v 1 the transferred energy can be approximated as 
AE = h(u — u') ~ hp, the transferred momentum then about Ap ~ hu/c. From 
this we can estimate the change of the kinetic energy as AE^ n ~ (Ap/rrii) = 
\h 2 v 2 1 (m;c 2 ). For a typical UV radiation \h 2 v 2 / (to;c 2 ) <C AE, and this 
means that most of transferred energy goes finally to heating (for v > v') or 
cooling (for v < u') of material and not to the macroscopic kinetic energy. This 
explains why the most common interaction of material and light results in heat- 
ing or cooling. On the contrary, in the case v ~ v' most of transferred energy 
goes to the macroscopic kinetic energy and the irradiated material is accelerated 
(assuming isotropic emission). 

Most effective processes for wind acceleration are those for which v ~ v' . 
This condition is fulfilled for light scattering in lines and on free electrons. Both 
processes are important for the acceleration of hot star winds. 

However, even line scattering may become less efficient. This is basically 
connected with the fact that frequencies of absorbed and emmited line radi- 
ation are slightl y different. This effect, introduced in the domain of hot star 
wind theory bv iGavlev &z Qwockil (^994) is called Doppler or Gayley-Owocki 
heating/cooling. An example of this effect is given in Fig. [21 Let us have an 
artificial light source that emits radiation only in a very narrow wavelength in- 
terval corresponding only to the left (blue) part of the line-profile (see Fig. |2J). 
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After the processes of absorption and emission the light from this source is redis- 
tributed over all wavelengths of a given line. Some part of the radiative energy 
has been thermalized, emitted radiation has lower energy, and the plasma can 
be heated by this process. On the other hand, let us have a light source that 
emits radiation with wavelengths corresponding only to the right (red) part of 
the line-profile. Again, the radiation is after the processes of absorption and 
emission redistributed over all wavelengths of the line. Clearly, after this, radi- 
ation has more energy that is taken from particle kinetic energy. Consequently, 
the plasma is cooler no w. T hese processes are impo rtant for low-density winds 
Csee lCavlev fc OwockilHoP : iKrticka fc KubatH200lT ). 




X X 

heating cooling 
Figure 2. Schematic picture of Doppler (Gayley-Owocki) heating/cooling. 



2.2. How to transfer momentum to the passive component? 

Since hot star winds are ionised, the most efficient way to transfer acquired 
momentum from heavier elements to the passive component (H, He) is due to 
the Coulomb collisions. Frictional force on passive component (p) due to metallic 
ions (i) is given by 

4Trq 2 q? 

/ P i ~ ripUi ^ 1 In AG(xip), (1) 

where n p , n\ are number densities of wind components, V{, v p are their radial ve- 
locities, and q p , qi their charges. The frictional force (J2) depends on the velocity 
difference v\ — v p via the so called Chandrasekhar function G(xi p ) (see Fig. |3J), 
where 

o? m ~2kT— ^. (2 



Q: 



For very low velocity differences, x; p < 0.1, the transfer of momentum between 
metallic and passive wind component is efficient. Wind is well-coupled in this 
case and it can be treated as one component. For hi gher velocity differences , 
Xi p 0.1, the frictional heating becomes important (|Krticka fc KubatJl200lh . 



For even higher velocity differences x- ip > 1, the Chandrasekhar function G(x- ip ) 
is a decreasing function of the velocity difference. Consequently, the transfer 
of momentum between metallic and passive wind components is inef ficient and 
wind components may decouple (e.g.lSpring mann &z Pauldrachlll992l ). 
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Figure 3. The Chandrasekhar function G(x- lp ). 



3. Passing to macro: the radiative force 



The radiative force is given by an integral 



/md = \ j X» F " dzA ( 3 ) 

Whereas the absorption coefficient \v depends explicitly only on local wind 
properties, the radiative flux F v is a non-local quantity and has to be obtained 
using the solution of the radiative transfer equation. Another complication for 
the calculation of the radiative force arises due to the Doppler effect. As the 
stellar wind is accelerated, the wavelength of a given line at which the wind 
absorbs the stellar radiation shifts to shorter wavelengths (in the reference frame 
of the star). Although this essentially enhances the radiative force since the wind 
is able to absorb stellar radiation that is unattenuated by the material below, this 
effect also significantly complicates consistent calculation of the radiative force. 
However, in the case of the stellar wind with relatively high velocity gradients 
both these effects may h elp to calculate the radiative force in an approximate 
way. This is the so called ISobolevI (j!947f > approximation . 

The radiative force due to optically thick line in the Sobolev approximation 
is given by 

jrad, thick ^ vF v 

c 2 dr ' 

where F v is the stellar flux at the frequency of a given line v and v is the wind 
velocity. Note that the radiative force in an optically thick line does not depend 
on line properties (e.g. the occupation numbers of corresponding levels or the 
line oscillator strength). This is one of the reasons why metals that have very 
low number density compared to hydrogen or helium (but a large number of 
optically thick lines) may be so important for the wind acceleration. 
The radiative force in an optically thin line is given by 

f ad, thin „ ™L Fvg . f ( 5 ) 

Metals dominate also to the optically thin line force due to very high number of 
their lines and due to frequent complete ionization of hydrogen and helium. 
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The total radiative force is given by the sum of contributions @ and © of 
individual lines. The calculation can be simplifie d using two approaches. First, 
it is p ossible to use the line distribution function ((Castor et al.lll975l : iPuls et all 
2000) parameterised by the set of force multipliers k, a and 8 to obtain the 
radiative acceleration in CAK approximation 



where p^ is electron density and dv/dr i s the ve l ocity gradient. Another line- 
force parameterisation was introduced bv iGavlevI (jl995T ) using Q parameter, 

rad nl-a J /ldiA a , . 

9 ~ Q Pcl {-pd^J ■ (?) 

4. Macro view 

4.1. Basic structure of ID CAK models 

We have understood how the stellar winds of hot stars are accelerated on a micro- 
level. Now we shall discuss the influence of the wind microscopic structure on 
the observable macroscopic properties. 

To do so, we have to calculate at least approximate stationary spherically 
symmetric wind models in the CAK approximation. However, th ere is a p r oblem 
because an extensive table of line-force parame ters calculated by [Abbott! (1982) 
leads to overestimation of mass loss rates (e.g. iKrticka fc Kubat!l2004al ) On the 
other hand, f or massive stars it is possible to use recent extended set of force 
parameters of iKudritzkil ( 2002f ). 



It is possible to show that wind mass loss rate scales with stellar luminosity 

M = Airp(r)v(r)r 2 ~ L 1 /"', a = a - 5. (8) 
( Kudritzki fc Puis! [2000), where a and 5 are usual CAK parameters. This 



means that wind mass loss rates depend mostly on the stellar luminosity. Mass 
loss rates are also significantly influenced by wind properties on a micro-level 
(e.g. metallicity, ionization structure) via the parameter a' . 

Wind terminal velocity Voo depends mostly on the escape velocity 

Voo = c(T ef[ )v esc , c(T cff ) « 1 - 3 (9) 

(e.g. lLamers et al.l[l995h and only slightly on wind properties on a micro-level. 

It can be also shown that clever combination of the mass loss rate, the 
terminal velocity and the stellar radius in the form of 

M Voo (R*/R G ) 1/2 , (10) 

that resembles wind momentum and is therefore called the modified wind mo- 
mentum depends mostly on the stellar l uminosity and only marginally on the 
stellar mass fe.g. lKudritzki fc PulslEiooT ) . 

The approximation described above allows calculation of wind models. De- 
rived velocity structure (e.g. Fig. |1J) can be in many cases approximated by the 
so-called beta velocity law using parameter (5 in the form v{r) = (1 — R*/r)^ . 
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Figure 4. Calculated velocity structure for 06V star in the CAK approxi- 
mation. The overplotted beta-velocity law (dashed line) is nearly identical. 



4.2. Multicomponent effects 

For low-density winds the transfer of momentum between metals and hydrogen 
or helium may become inefficient (see Sect. I2.2T1). There are several ty pes of flow 
with respect to the multicomponent effects ((Krticka fe KubatJ|2004bl h 

Winds with negligible multicomponent effects For stars with relatively dense 
winds (e.g. winds of Galactic O supergiants) the multicomponent effects can be 
neglected. Such stellar wind can be adequately treated as one-component. 



Wind temperature influenced by frictional heating For stars with lower-density 
winds or very low metallicity the transfer of momentum (and energy) between 
metals and passive (H, He) component becomes ine fficient and part of transferred 
energy goes to heating (|Cure1ll992l : IKrticka et alJl2003th 

Decoupling in the wind For stars with very low wind densities or with very low 

metallicities hydrogen and he lium decoupling may occur in the wind (ISpringmann <fe Pauldrachl 
1992; Krticka &: Kubat 2001). The stellar wind is not stable in this case (lOwocki fc Pulsl 
200 21: IKrticka fc Kubadl2002l). This problem was studied using HD simulations 
flPorter fc Skouzalll999l : IVotruba et alll200fil ) . 

Decoupling of wind c omponents in the atmosp here Helium decoupling in the 
wind was proposed bv lHunger Ik Grootd ( 19991 ) as the explanation of the chem- 



ical peculiarity of He-strong stars. Due to its low charge, helium may decouple 
in the stellar atmosphere of cooler stars and a helium-free wind may exist in 
this case. For extremely low-density winds also hydroge n may decou ple in the 
atmosphere and purely metallic stellar wind may exist (Babelllfl996). For hot 



stars with lower luminosities the radiative force is not able to expel atoms out 
from stellar gravitational potential well, however the radiative force may cause 
chemical peculiarity (e.g. lMichaudll2004l ) . 



The regions in HR diagram with different types of stellar wind are given in 
Fig. |SJ Note, however, that NLTE models are necessary to study these effects in 
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Figure 5. Regions with different types of stellar wind in HR diagram. 



detail. Finally, multicomponent effects may become important even for slightly 
earlier stars due to possible decoupling of individual elements. 



5. Beyond classical CAK models 
5.1. Influence of rotation 

Rotation may cause significant deviation s of wind structure from the spherical 
symmetry. iBiorkman fc Cassinellil (^93) proposed that discs of rapidly rotat- 
ing stars may be caused by the wind com pression due to the stellar rotation. 
Howe ver, 2D hydrodynamic wind models ((Qwocki et al.l Il996l : iPetrenz Sz Pulsl 
2000) showed that due to the gravity darkening and nonnegligible non-radial 
component of the radiative force the wind density and mass loss rate at the 
equator are lower than at the poles. 



5.2. Influence of magnetic fields 

Also magnetic fields may cause significant deviations from spherical symmetry. 
They may be i mportant for the correct explanation of the circumstellar activ- 
ity of Bp stars (|Babel fc Montmerldll997l: iGroote Hungerlll997l : iTrigilio et all 
12004ft . According to MHD models of lud-Doula fc Owockil fe002ft . the overall de- 
gree to which the wind is influenced by the magnetic field depends on the ratio 
between magnetic field energy density and wind kinetic energy density 77*, 

pv z /2 

For a weak confinement (77* < 1), the magnetic field is opened by the wind out- 
flow. The structure of the circumstellar magnetic field is given mostly by the 
stellar wind and the influence of the magnetic field on wind is not a significant 
one. However, for strong confinement (77* > 1) the situation is essentially oppo- 
site. The flow near the star is driven by the magnetic field. For B stars even a 
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moderate magnetic field intensity (B < 100 G) causes strong confinement of the 
circumstellar environment by the magnetic field and is therefore very important 
for its structure. Note that for very stron g magnetic fields some par ts of the 
circumstellar envelope may be hy drostatic (|Tow nsemj_fc Qwockill2005t ). 

Magnetic fields were used bv lCassinelli et al, (2002) to ex plain the Be phe - 
nomenon. However, this is not supported by MHD simulations (|ud-Doulal l2006 ) . 

5.3. NLTE models 

The "triplet" of parameters (k, a, 5) or (Q, a, 5) gives only rough approximation 
to the radiative force. To solve this problem, it is possible to either 

• introduce depth-depe ndent radiative- force parameters fe.g. lPauldrach et alJ 
120011 : lKudritzkill2002ft . or to 

• calculate the radiative force directly without any radiative-force parame- 
ters fe.g. lVink et al.l200"ll : lKrticka fc Kubatl2004al . and references therein). 

NLTE approach is necessary in any case to obtain correct wind parameters 
(mainly mass loss rates). 

The introduction of a more realistic radiative force based on the appropriate 
level occupation number s enables detailed study of wind parameter variations. 
iPauldrach h Pulsl (|1990T ) found high sensitivity of cal culated wind para meters 
of P Cyg on its stellar parameters - the hi- stability. IVink et alJ ( 1999h found 



the bi-stability jump at around T P g 25 000 K for no rmal supergiants. Using 
terminal velocity measurements of lLamers et al.1 (|l995h they concluded that for 
stars cooler than the temperature corresponding to the bi-stability jump the 
mass loss rate M increases 5 x , whereas the terminal velocity Voo decreases 2 x . 
The bi-stability jump is caused by an increase of the line acceleration due to 
Fe III lines close to the stellar surface. These calculations slightly overestimate 
the correct temperature location of the bi- stability jump , which occurs roughly at 
T e fj « 21 000 K according to the results of lLamers et al.1 ( 1995T). Jump properties 
have s till to be tested against observations since results of [Trundle Sz Lennqnl 
(2005) show a different picture of the jump for SMC stars (c.f. lSearle et alJl20Qq ). 

5.4. Radiative force — more exact approximations 

Inclusion of higher order app roximations to the radiative force leads to the wind 
instabilities (Fel dmeierlll998l ). Resulting wind shocks may (at least partly) ex- 
plain the observed X-ray emission o f hot stars. Det ailed discussion of this prob- 
lem can be found e.g. in a review of lOwockil (j2004h . 



6. Open questions 

There are many open questions connected with the stellar wind of OB stars. 
Here we discuss only those that seem to as to be especially appealing nowadays. 

6.1. Influence of instabilities and clumping 

Observed hot star wind properties show signatures of spatial l y org anised struc- 
tures - clumping (e.g. lAntokhin et al.1 fl988l : iMartins et al.1 120051 ) . As noted 
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above, the ra diative driving is unstable and may cause generation of shocks (e.g. 
IOwockill2004l ). The relation between instabilities and clumping is still unclear, 
although clumping parameter is ofted used as an additional free parameter of 
wind models. 

6.2. Winds close to limit 



Durin g the stellar evolution some stars may come close to the limit ( Maeder Mevnet I 



2000). However, it is not clear what sets the mass loss rate in this case. The 
radiatively-driven wind cannot drive outflow with an arbitrary mass loss rate 



( Owocki &: Gavlevlll997f ). It is possible that to estimate these mass loss rates it 



is necessary to calculate "unified" models of stellar interior and envelope. 
6.3. What are wind mass loss rates? 

Reliable theoretical predictions of the mass loss rate are still insufficient. One 
of the reasons is that detailed (NLTE) wind models are necessary to predict 
mass loss rates. As the consequence, reliable predictions (for OB star domain) 
are available only for O-stars and luminou s B-stars (e.g. IPauldrach et al.ll200ll : 
IVink et al. Il200ll: iKrticka fc Kubatl l2004aT ) and for hot horizontal branch stars 
(IVink fc Cassisill2002D . Note that even frequently used predictions of Vink et al 



200lh suffer from many simplifications like neglect of "line branching" rtSim 



2004), wind instabilities and clumping or X-ray radiation ((MacFarlane et al 



1994). Whereas there is a relatively good agreement between theoretically pre- 



dicted mass loss rates and mass loss rates derived from observations for indi- 
vidual hot OB stars, the re is a significant disagreement between these values 
for cooler B supergiants ( Vink et al.ll2006T ). Although part of this discrepancy 



may be due to observations, some part of this discrepancy is likely due to model 
simplifications. 

Even worse, for many stars (e.g. for many main-sequence B stars) there are 
no reliable predictions of mass loss rate available. The problem is that there are 
several processes that may influence the mass loss rate of these stars, e.g. 

• multi compo nent wind structu re dSpringmann Sz Pauldrachll992l : ICure1l992l : 
lBabellll995l : iKrticka fc KubatlfcOQlh . 

• GO (Doppler) heating dGavlev fc Owockil 1 1994ft . 



thin- wind effect (lOwocki fc Pulsl Fl999 N ). 



Low-density stellar winds are also very difficult to observe. However, some 
indir ect wind indications may be available, e.g. due to c ertain chemical peculiar- 
ities ( Landstreet et al.lll99S:lDworetskv BudafcOOOT) or due to wind magn etic 



braking of stellar rotation ( Mikulasek et al.ll200 A lOksala fc Townsendl EoOBV 



There is an incre asing observational evidence that hot star winds are clumped. 
iMartins et all ((2005) studied winds of several Galactic O-stars in detail and con- 



cluded that mass loss rates derived assuming clumped winds may be 2 x — 5x 
lower than those derived assuming smooth winds (Fig. EJ). Finally, for some 
stars with low luminosities the wind mass loss rates seem to be more than ten 
times lower t han the predicted o nes (Fig. EJ). Simil ar effect was detecte d for 
SMC stars bv lBouret etaTI ((2003). Thin- wind effect (lOwocki fc Pultfi fafn may 



help to understand the origin of this difference. 
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M[M yr~ 1 ] calculated - Vink et al. (2001) 



Figure 6. Comparison of mass loss rate derived from observations allowing 
for wind clumping and the theoretical (smooth) mass loss rates 



7. Conclusions 

We have discussed basic physical principles that drive the stellar winds of OB 
stars both on micro and on macro levels. We have also seen that despite the 
tremendous effort of many astronomers, more than 30 years after the publication 
of the seminal CAK paper we are still not sure what the mass loss rates of hot 
stars are. To conclude wi th something more optimist ic we switch to slightly dif- 
ferent stars - WR stars. iGrafener &z Hamannl ( 2005h were able to theoretically 



explain the acceleration of the stellar wind of WR stars. According to these mod- 
els, wind sonic point is located deep in the stellar interior and the stellar wind in 
the inner regions is accelerated due to the iron opacity bump and convection in 
WR stars occurs also as a consequence of the iron o pacity bump. Iron opacity 
bump stands also beh ind pulsations in (3 Cep stars ((Moskalik Dziembowskil 
Il992l : ICox et all I1992T ) . Consequently, processes of the same physical nature 
cause completely different behaviour in different circumstances! Since wind crit- 
ical point in the case of WR stars is located relatively deep in the stellar interior, 
any differentiation between stellar core, atmosphere, and wind may become ar- 
tificial and we might need to (at least for some specific problems) calculate 
"unified" models of stellar interior and exterior. 
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Owocki: I must disagree with the notion that Fe opacity bump can be an 
important mechanism for wind driving. I think to focus on using this bump to 
set the wind through the sonic point is mis-guided. The sound speed is of order 
25kms~ 1 , about 30 times smaller than the escape speed v esc (~ 700 km s -1 ); 
implying an energy that is 1/1000 of that needed to escape the star Thus even if 
Fe opacity sets the wind through the sonic point, we must appeal to some other 
mechanism to do the 99.9% of the remaining work needed to escape. I think 
understanding wind mass loss should focus instead on the mechanism that does 
this overwhelming bulk of the work. 

Krticka: I agree that iron bump opacity alone is not able to drive the stellar 
wind of WR star and does only a certain part of work necessary to launch the 
wind. Also other opacity sources should be included. To drive the wind it is 
necessary to accelerate the wind material from the subsonic velocities to the 
velocity higher than the escape speed. 

Townsend: Another way of thinking about the iron opacity bump is to recall 
that it drives SPB & (5 Cep pulsations. The reason these stars show pulsations is 
that the iron opacity bump disappears when the temperature varies significantly 
from 200 000 K. This makes it very difficult to understand how a wind - which 
requires a spatially-extended region of high opacity - could be driven by the iron 
opacity b ump. 

Krticka: iGrafener Sz Hamannl ((2005) were able to consistently accelerate the 
wind from the subsonic velocities to the velocities higher than the escape speed. 
I do not think that your argument is important, because the physical mechanism 
that drives the pulsations is different. Pulsations are driven basically by the heat 
flux, however the stellar wind is driven by the gradient of the radiation pressure. 
Stee: The terminal velocity and the mass flux is strongly depending on the 
inclination angle (on the stellar colatitude) i.e. larger at the poles and lower at 
the equator. Thus it is normal that there is a large scattering in the mass loss 
rate versus effective temperature graph. Moreover the mass loss also depends on 
the lines you are using to determine it (IR, visible, UV). Finally, it seems difficult 
for me to compare observational and theoretical Voo and M for nonspherical wind 
(for instance following the bi-stability scheme you have shown). 
Krticka: Rotation is important for wind structure, however it has only sec- 
qnd order effect for st ars with rotational velocities well bellow the critical one 
( Petrenz & Pulsll200Cl L It is not likely that it causes order of magnitude differ- 
ences of M for cool B stars or too high predicted M for low-luminosity stars. 



